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ABSTRACT 


Wind-off,  bench  tests  on  a  wind  tunnel  model  that  had  freedom  to  pitch, 
and  yaw,  but  no  freedom  to  roll  indicated  that  the  pitching  and  yawing  motions 
were  inertially  coupled.  An  order  of  magnitude  analysis  showed  that  this 
coupling  was  due  to  the  product  of  inertia  term  The  inertially  coupled 

equations  of  motion  are  solved,  and  a  method  is  presented  for  determining 
the  magnitude  and  sign  of  from  recordings  of  the  model  motion.  The 

method  used  to  dynamically  balance  the  model  is  described,  and  experimental 
results  are  presented  that  verify  this  method. 
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NOMENCLATURE 


^e’  \ 

I  ,I  .1 

x’  y  z 

I  ,'I  ,I 

xy  xz  yz 
M  ,  M  ,  M 

X  y  z 

p.  q.  r 

qi  .qz 
t 

X,  y,  z 
x^, y^, z^ 

xjy,  0,0 

^0’ V 
<^11.012.0 

^6’ ^4/ 


Amplitude  components 
Flexure  damping  constants 
Constant  defined  by  Eq.  (19) 

Moments  of  inertia  about  body  axes 
Products  of  inertia  about  body  axes 
Moments  about  body  axes 

Angular  velocity  components  on  body  axes  system 
Normal  coordinates,  Eq.  (A-27)  and  (A-28) 

Time 

Body  axes  system  (Fig.  1) 

Space  fixed  axes 

Constants  defined  by  Eqs.  (A-38)  and  (A- 39) 
Eulerian  angles 

Constants  defined  by  Eqs.  (14)  and  (15) 

Constants  defined  by  Eqs.  (16),  (17)  and  (18) 
Circular  frequency  in  pitch  and  yaw 


I.  INTRODUCTION 

In  order  to  investigate  possible  aerodynamic  interactions  between 
the  pitching  and  yawing  motions  of  a  nonrolling  body  undergoing  non- 
planar  motions,  a  flexure  was  designed  to  give  a  wind  tuniiel  model 
freedom  to  either  pitch,  yaw  or  pitch  and  yaw  simultaneously.  The  model, 
a  10°  half  angle  cone,  was  weighted  so  that  the  body  axes  of  rotation 
were  intended  to  be  the  principal  axes  of  inertia,  and  yet  ly  was  not  equal 
to  1^.  Releasing  a  model  balanced  in  this  manner  from  a  combined  pitch 
and  yaw  attitude  with  zero  initial  velocities  results  in  a  nonplanar 
Lissajous  motion,  that  is,  a  motion  with  different  frequencies  in  pitch 
and  yaw.  As  the  model  is  not  rolling,  the  pitch  and  yaw  components  of  the 
motion  should  be  uncoupled,  damped,  sinusoidal  oscillations. 

Wind-off  bench  tests  with  the  model  in  vacuum  did  not  give  the  uncoupled 
motion  that  was  expected.  The  pitching  and  yawing  motion  envelopes  were 
not  smoothly  damped  exponentials  but  rather  "  sinusoidally"  damped  expo¬ 
nentials.  The  sketches  below  are  an  exaggeration  of  the  motion  obtained. 

•  > 

The  initial  velocities  in  both  sets  of  sketches  are  zero  (  ^  =  0),  with 

Ip  Q  =  0Q  in  the  first  set,  and  4^  q  ~  0,  in  the  second  set.  This  beating 
type  of  motion  indicates  that  the  pitch  and  yaw  motions  are  coupled  due  to  a 
slight  inertia  asymmetry,  that  is,  the  principal  axes  of  inertia  do  not  coincide 
with  the  axes  of  rotation.  The  purpose  of  this  report  is  to  determine  which 
product  of  inertia  causes  the  coupling,  to  determine  its  magnitude  from  the 
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II.  THEORETICAL  ANALYSIS 


The  moment  equations  of  motion  for  a  body  fixed  coordinate  system 

(Fig.  1),  with  its  origin  located  at  the  c.  g.  of  the  body,  can  be  written  as 

M  =  I  p  +  (I  -  I  )rq  -  I  (r  +  qp)  +  I  (rp  -  q) 

X  X  z  y  xz  “  xy 


(1) 


+  I  (r^-  q") 

yr 


(1) 


M  =  I  d  +  (I  -  I  )pr  -  I  (p  +  qr)  +  I  (pq  -  f  ) 
y  y  X  z  xy  yz 


(2) 


M  =  I  r  +  (I  -  I  )qp  -  I  (q  +  pr)  +  I  (qr  -  p) 
z  z  '  y  X  yz  ^  ^  x:z  ^  ^ 


+  I  (q2  -  p') 
xy  ^  ^ 


where  p,  q,  r  are  the  components  of  the  angular  velocity  along  the  body 
axes,  X,  y,  z  respectively.  The  angular  velocity  components  can  be  written 


(3) 


as 


ps  (j)  -  x(y  sin  0 
q  =  0  cos  (j)  +  xf/  cos  9  sin  6 

r  =  ij/  cos  9  cos  (j)  -  9  sin  (j) 


(4a) 

(4b) 

(4c) 


The  terms  M^,  M^,  M  are  the  moments  acting  on  the  body.  For  a  nonrolling 
vehicle  in  a  vacuum,  without  wind,  these  terms  represent  the  flexure  moments 
acting  on  the  system.  The  pitch  and  yaw  moments  can  be  written  as 
M  = 


C^0  -  k^0 


^z  = 


(5a) 

(5b) 
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For  0:  =  0=0  and  assuming  small  angles,  Eq.  (4)  reduces  to 
p  Ri  -  xjy  9 
q  .=  6 

'  r  ;  W  0/ 

Substituting  Eq.  (6)  into  Eqs.  (-2)  and  (3)  gives: 

;,M  :  =  :I0+(I  -I  )(-0/0^)-  I  {-xi/e)+l  (-0/0  0-0^) 

'  Y  '  5c  z  3C  y  2 


-  +  .I  (-0^02-0/) 


xz 


bM„.=  ;I„0+(I  -I  )(-0/00)-I  (0-0/0/0) 


y  X 


yz 


(6a) 

(6b) 

{6c) 


(7a) 


.;+!  (0  0/+  0/0  +  0/0)  +  I  0^-x{/e^) 

■  xz  xy 

■jT.bgset equations  can  be  simplified  due  to  the  small  angle  assumption.  Assume 

i/that  p.itch  £  and -yaw,  angles  can  be  approximated  as 

(9  ~~  ^9^  cos  Wgt,  0/  =  0/^  cos  cu^t 


(7b) 


o 


iTh.en, 


.and 


!  0r-  =  .  -  0  sin  C0_t  ,  0/  =  -  O),  0/  sin  W,  t 

'  0  o  0  ^  0/^0  0/ 


^  0r  =-  -  d>50  cos  CO  „t  ,  0/  =  -  (of  0/  cos  C0,t 

0o  0  0/’^o  0/ 


AArt  order  of  niagnitude  analysis  can  now  be  made  for  the  individual  terms  of 
EEq.  17) 
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i^/  0  w  li/  0  0)^ 

^^oo\p 

•  • 

10  0  fi  10  0  a),a)_ 

^  ^o  o  \y  6 

6'^  6^  "dj  9^00,  CO. 

^  ^o  o  9 

,  ;0^  0  «  i//^0  co^ 

^  o  O  \f/ 

102  S»  10^00^, 

^  ^o  \l/ 

9^  «  0^co^„ 

o  0 

0  «  0  co^ 

o  0 

xU  ^  \l/o/, 

^  ^o  Tp 

d/^  6^  ~  p^9^oj\ 

^  ^o  o  p 

•  •  •  * 

All  the  product  terms  are  small  compared  to  either  the  9  or  p  term  if 


the  amplitudes  p  and  0^  are  small.  Thus  for  small  angles,  Eqs.  (7a) 
and  (7b)  reduce  to 


M  =  I  0’  -  I  p’ 

y  y  yz 


(8a) 


M  =  1  p  -  1  9 

z  z  yz 


(8b) 


As  can  be  seen  from  Eqs.  (8a)  and  (8b),  the  product  of  inertia  term  I 

yz 

causes  the  coupling  between  the  pitch  and  yaw  motions.  Although  I  and 

yaz 

I  may  be  of  the  same  order  of  magnitude  as  I  ,  they  are  multiplied  by 
xy  °  yz 

relatively  small  quantities  so  that  they  can  be  ignored  compared  to  the  terms 
I  0  and  Ip. 

yz.  yz 

The  effect  of  I  on  the  pitch  and  yaw  motion  and  the  value  of  I 

yz  yz 

causing  the  coupling  can  be  determined  from  the  solution  of  Eqs.  (8a)  and 
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(8b).  Substituting  Eq.  (5)  into  Eq.  (8)  gives 


I  0-1  =  -  C.  0-  k_  0 

y  yz^  0  0 


I  1//-I  e  =  -  C,  -  k  , 

z  ^  yz  l//^ 


(9a) 

(9b) 


oi',  in  matrix  form, 


I  -I 
y  yz 

-I  I 
yz  z 


^0  ° 

°  % 


%  ° 
«  V 


(10) 


Equation  (10)  is  a  second  order,  coupled,  differential  equation.  As  the 

coefficient  matrices  are  symmetric,  this  equation  can  be  solved  using  a 

(2  3) 

normal  coordinate  system'  ’  The  initial  conditions  used  to  solve  Eq.  (10) 
are 


=  0 


=  0 
o 


•j^(0)  * 


(11) 


Ediio'wing  a  method  identical  to  that  given  in  Ref.  3,  the  solution  to  Eq.  (Ip) 
is  derived  in  A-ppendix  A .  The  results  are 


-Ant  )  0 


Iz 


(0,1  '/'o  4))  ,  ^0  , ,  , 

— -  +  —  (0  xl/  -  9) 


w 


0 


0 


d-  (B  |■%“  (%"*/’  l// )cosX0,t 

i(  0  °  TzTi'  ,0/  0 


(4-0,  >o)  ^0/ 

o),  w,  '  o  q2T>' 

0/  0/ 


sin  Cdgt 


sinco,  t 

0^ 


(12) 


and, 


n&sco,  t  -+  ~ 

.0/  0 


- L2 -  ^  (0  _0  0/  ) 

O),  W,  o  tzh) 

0/  0/ 


■"'^4  1  (0,, ,  ,  1 

G  ]  /-w  COS  CO^t 

(9  0  0 


(0n’4-4)  .  ^  _ 

-  +  —  (0  xb  -  9  ) 

CO.  co^  '  ii  o' 


0 


0 


sinco,  t 

0^ 


sino^tj 


(13) 
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where 


= 


X,  = 


0  = 

hi 


(!)  = 

hz 


V.h "  ^ 

<‘A‘'0  ^  V 


C  o)^ 

0  0 


C^  +  C  ,  (J,  Q,  a)f. 

VJ _ !L  Jl  ~  jtjk 

0^  k^  +  k  ,  2  2k  , 

r.  0  \p  i/y 


I  CO" 
vz  I 


I  CO"  -  k 
z  tj 

T 


I  CO"  -  k  I  CO  , 

y  ijy  0  yz  xl/ 


I  CO"  I  co^  -k 

yz  0  _  z  0  ! 

I  -  ^3  ""  I  CO^ 
y  0  0  yz  0 


0-0 

Til  h 


= 


I  I  -  (I  )' 
y  z  yz 


k^  I  +  k  ,  I 
0  z  0/  c 

21^ 


.  k.I  +  k,  I 

10  2;  0/  V 


k.  k, 
0  0/ 


The  sign  in  front  of  the  square  root  term  in  Eq.  (20)  is  determined  from  the 

case  of  zero  I  ,  where  (of  =  k,  / 1  ,  and  CO^  =  k^/l  .  If  k,  /l  is  greater 
yz  0/0/  z  00y  0/z‘^ 

than  then  the  plus  sign  in  Eq.  (20)  will  be  associated  with  to^  ,  and 

the  minus  sign  with  CO^.  If  k  /I  is  greater  than  k,  /I  ,  then  the  reverse 


0'  y 


will  be  true. 


The  effect  of  I  on  the  model  motion  can  be  seen  from  Eqs.  (12), 

yz 

and  (13).  As  there  are  two, nearly  equal,  frequency  components  present  in 
each  equation,  the  resulting  motion  for  both  pitch  and  yaw  will  be  a  beating 
type  of  motion.  That  is,  the  envelopes  will  have  a  sinusoidal  type  pattern. 


7 


<;%. » ....:i.mr.,->i'L-.^-j^o^vg.^t‘l..  r=aa.^^-,  r.- 


,The  effect  of  I  on  the  frequencies  can  be  determined  from  Eq.  (20), 
which  can  be  rewritten  as 


(%)■■ 


(I  /I  +  k,  /Ic  )  1/1  +k,/k.  2 

z  y  xj/  6  ^  z  V  W  6 

2(1  /I  -  /l^)'l.  2(1  /I  -I^  /i2 

z  y  yz  y  z  y  yz  y 


k  I.  /k^ 


(I  /I  -  /I^) 

z  y  yz  y 


”2(1  /I  -  /I^) 

z  y  yz  y 


/I  + 


(I  /i  -1^/1^ 
z  y  yz  z 


(I  /I  -  I  ^  /i^) 
z  y  yz  y 


Equations  (21)  and  (22)  are  plotted  in  Figs.  2  and  3,  as  the  frequency  ratios 

tO«/(k^/I  )v^  and  O),  /(k  /I  )'/2  versus  I  /l  for  k  /k  equal  to  1.014, 

'  ,0/ i.' 0/ :  y  ai.  .  .  i//  '  x(/'  z  yz  y  x/y  6 

and  for  several  values  of  I  /  I  ,  These  curves  show  that  I  has  a  small 

■'  e  5  f.;-  £  er  c.  V :  .  -  ■  z  y  ,  yz 

Ojffect  on  the  frequency,  causing  it  to  increase  or;decreasc  depending  on  the 

value  of  I  / 1  . 

'  ■  -  c:  iZ,  _y. 

The  moments  of  inertia,  I  and  I  ,  as  well  as  the  product  of  inertia, 

■j  !.e  ir.  y  z 

I  ,  can  be  obtained  from  Eqs.  (16)  to  (18)  as 

■  yZ:  c  e  J-;  L  :  u  .  -  : 

k  (j)  ool  -  (j)  oo] 

0  hi  0  12  d/ 


1  ~ 

-Z  " 


-S:'  /  ,2  /  ,2 
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In  order  to  determine  the  I  of  a  system,  values  of  k,  ,  6  ,  00 , 

yz  ^  xj/  ^  if/ 


and  OOq  are  required.  While  the  flexure  stiffness,  known  from 

static  bench  tests,  ^q>  ^  must  be  determined  from  dynamic  tests, 

i.  e.  ,  from  recordings  of  the  pitching  and  yawing  motions. 

Values  for  0^^  and  0^^  can  be  determined  as  follows.  Assuming  that 
the  initial  angular  velocities  are  zero,  Eqs.  (12)  and  (13)  reduce  to 

-  Aflt 

Wi.'t.-  "fl*  + 

e 


and. 


0/(t)  =  e  {9^-  0^,0/^)sinCO^t 


-V  (0n'^6-0o) 


+  e  ^ 


^9 


- cosco^t  +  ^ - (0,,  \p  -  9  )sina^„t 

0  9  0a).  '^11^0  o'  9 


or. 


where 


9{t)  = 

-  A^t 

A.  e  cos(a)  t 

9  9 

-  A0,t 

(26) 

0/(t)  = 

-X  ,  t 
ih 

A,  e  cos(a),  t 

0/  0/ 

-  A^t 

(27) 

^9 

12  0 

« 

fA 

^9 

(28a) 

^  <«0-  V 

V 

(28b) 
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therefore,  from  Eqs.  (Z8)  and  (29) 


(29a) 

(29b) 


h2  B, 


The  sign  of  I  is  associated  with  the  sign  of  (j)  which  depends  upon 
yz 

the  signs  of  the  A's  and  B's.  If  the  envelope  of  the  motion  initially  decreases 
then  A  and  B  in  Eq.  (26)  or  (27)  will  have  the  same  sign.  However,  if  tiie 
envelope  initially  increases,  then  A  and  B  will  have  opposite  signs.  The 
sign  of  A  always  has  the  same  sign  as  the  initial  release  angle. 


(30) 


III.  EXPERIMENTAL  INVESTIGATION 

The  objective  of  the  experimental  investigation  was  to  determine  the 

-I^^  of  the  model,  and  then  to  eliminate  it  by  properly  balancing  the  model. 

The  I  of  the  original  model  was  small,  and  therefore  it  was  difficult  to 
yz 

determine  the  amplitude  components  of  the  motion.  In  order  to  overcome 

this  difficulty,  a  larger,  known  I  was  added  to  the  model.  The  resulting 

yz 

motion  was  more  easily  reducible.  The  original  I  of  the  model,  plus  the 

yz 

added  I  was  determined  by  fitting  the  data  to  Eqs.  (26)  and  (27)  by  the 
yz 

(4  5) 

least  squares  differential  corrections  technique'  ’  .  The  added  I  was 

then  subtracted  to  give  the  I  of  the  model.  This  was  repeated  with 

yz 

another  known  I  ,  and  the  model  I  was  taken  as  the  average  of  the  two 
yz  yz 

fpsults.  •  . 

PROCEDURE 

The  two  model  configurations  tested  were  the  basic  model  plus: 


Configuration 

I  sl-ft^ 

yz  added 

I  sl"ft^ 

V  added 

I  . ,  ,  sl-ft 

z  added 

(1) 

-5.  25  X  10-5 

4.  843  X  lO-'^ 

5. 393  X  10-4 

(2) 

1-4.  27  X  10-5 

4, 843  X  10"^ 

5. 393  X  10-4 

por  each  model  configuration  the  model  was  first  balanced  in  the  xz,_and 
yz  planes  to  insure  that  the  center  of  gravity  was  located  at  the  pivot  point. 
The  model  was  mounted  in  a  test  cabin  which  was  evacuated  to  less  than 
2mm  of  mercury  to  eliminate  aerodynamic  effects.  The  model  was  set  at  a 
predetermined  pitch  and  yaw  angle,  and  released. 
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The  resulting  angular  motion  was  recorded  by  strain  gauges  mounted 

on  the  flexure.  The  analogue  output  from  the  gauges  was  digitized,  and 

stored  on  magnetic  tapes.  Calibration  runs  for  the  conversion  of  the  digitized 

data  to  degrees  were  made  before  and  after  each  run.  The  pitch  data  was 

fitted  to  Eq.  (26),  and  the  yaw  data  to  Equation  (27).  The  section  length  fitted 

was  one  beat  period.  The  amplitude  components  determined  from  the  least 

squares  fit  were  then  used  to  calculate  Eq  (30),  and  then  0, 

Equation  (18).  The  value  of  <p,  along  with  Wq  ,  and  CO^,  which  were  also 

obtained  from  the  least  squares  fit,  were  used  to  calculate  I  ,  Eq.  (25), 

yz 

as  well  as  I  ,  Eq.  (23),  and  I  ,  Equation  (24).  Three  consecutive  sections  were 

y  z 

fitted,  and  the  average  values  resulting  from  these  three  fits  were  taken  as 

the  I  ,  I  ,  and  I  of  the  configuration.  The  added  values  of  I  ,  I  and  I 
yz  y  z  yz  y  z 

were  subtracted  from  the  measured  values  to  give  I  ,  I  and  I  of  the 

yz  y  z 

original  model. 

RESULTS 

The  motions  that  were  obtained  from  the  runs  are  shown  in  Figs.  4, 

5,  and  6.  As  can  be  seen  from  run  #1,  Figs.  4a  and  4b,  the  9  envelope 

initially  increases  while  the  l//  envelope  initially  decreases.  The  I  ^  of 

this  configuration  should  be  negative.  This  is  in  accord  with  the 

I  =  -5.  25  X  10  ^  sl-ft^,  that  was  added  to  the  model.  For  run  #2,  Figs.  5a, 
yz 

and  5b,  the  situation  is  reversed.  The  6  envelope  initially  decreases. 


while  the  xjy  envelope  initially  increases.  The  of  this  configuration 
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should  be  positive.  This  is  again  in  agreement  with  the  I  =  +4.  27  x  10"^ 

yz 

sl-ft^  that  was  added  to  the  model  for  this  run. 

The  results  of  the  least  squares  fitting  procedure  are  summarized  in 
Table  I. 


TABLE  I: 

Results  From  Least’ Squares  Fit  to 

Data 

Configuration 

%model 

^zmodel 

V^model 

1 

0. 002732 

0.  002397 

0. 30  X  10"^ 

2 

0. 002737 

0. 002405 

0. 37  X  10-5 

The  average  value  for  the  I  of  the  model  is  0.  34  x  10"^  sl-ft^  with  a 

yz 

probable  error  of  ±0.  11  x  10“^  sl-ft^,  or  about  32%  of  the  average  value. 
The  average  values  for  and  are  0.  002735  sl-ft^,  and  0.  002401  sl-ft^, 
respectively.  The  probable  error  of  these  terms  is  estimated  to  be  about 
0.  1%. 


The  angle  between  the  principal  axes,  and  the  body  axes  is  given  by 


1  -a  /2I 

r 


,  which  for  small  angles  reduced  to  --  ,  For  runs  1 

•^z”-^y 


and  2,  with  intentionally  added  I  ’s  (Figs.  4  and  5),  this  angle  is 

yz 

10.  2  degrees  and  9.  6  degrees  respectively.  For  the  original  model,  the 
angle  is  0.  6  degrees.  Even  this  small  angle  of  less  than  one  degree  was 
sufficient  to  cause  inertial  coupling  of  the  pitching  and  yawing  motion  that 


could  be  detected  on  oscillograph  records. 
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IV.  DYNAMIC  BALANCING  METHOD 


5  ? 

The  I  of  the  model  was  determined  to  be  0.  34  x  10"  sl-ft  . 
yz 

An  I  of  -0.  32  X  10'^  sl-ft^  was  then  added  to  the  model  by  placing  equal 
yz 

masses  at  45°  and  225°  on  the  model  (see  sketch).  This  should  reduce  the  net  I 


yz 


of  the  model  to  approximately  zero.  The  model  was  statically  balanced  by 

placing  masses  on  the  geometric,  or  principal  axes,  i.  e.  ,  0°,  90°,  180°, 

or  270°.  This  prevented  any  additional  I  being  added  to  the  model.  This 

yz 

configuration  was  run,  and  the  resulting  motion.  Figs.  6a  and  6b,  showed 
that  the  model  is  in  near  perfect  trim. 
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V.  CONCLUSION 


The  present  study  of  the  motion  in  vacuum  of  a  nonrolling  body  with 
freedom  to  pitch  and  yaw  shows  that  the  product  of  inertia  is  the 

dominant  term  causing  inertial  coupling  of  the  pitching  and  yawing  motion. 
The  presented  analysis,  and  the  described  method  of  eliminating  the  inertial 
coupling  have  been  verified  by  experiments.  The  model  is  now  ready  for 
wind  tunnel  tests,  the  results  of  which  will  be  reported  later. 
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VII.  APPENDIX  A 


The  governing  equations  of  motion  in  a  body  fixed  coordinate  systern 
with  its  origin  located  at  the  c.  g.  are 


I  0  -  I  rj/  +  CJ  +  kj  =  0 

y  yz  6  6 

-I  e  +  1  \}y  +C,j//  +  k,if/  =  0 


(A-1) 

(A-2) 


Consider  first  the  undamped  equations  of  motion 
I  0  -  I  4/  +  k^  6  =  0 

y  yz  0 

-  I  0  +  I  li/  +  k,;//  =  0 

yz  z  \}/ 

In  matrix  form,  Eqs.  (A -3)  and  ;{-A-4),are  written  as 


(A -3) 
(A-4) 


I  -I 
y  yz 

-I  I 
I  2 . 


+ 


0  k 


r0i  ^o■ 


[ijy)  »o. 


One  possible  solution  has  the  form 


0 


A2 


sin  OJt 


(■A-5) 


(A-6) 


Substituting  Eq.  (A-6)  into  Eq.  (A-5)  gives 


I  -I 
y  yz 

-I  I 
L  yz  z  J 


-rAi  or  s  incut 
-Azo/  sincut 


+ 


‘'9  “ 


0  k 


A I  s  incut 
A  2  s incut 


or, 


-I  CU^  I  OU^' 
y  yz 

I  CU^  -I  CU^ 

L  yz  Z 


■ 

■^1 

0 

_ I 

0 

+ 

_ 

■^z 

0 

Az 

0 

23 


matrix  addition  then  gives 


'k  -  I  0)2 

9  y 

I  0)2 

yz 

Ai 

0 

I  0)2 

L  yz 

k,  -  I  0)2 

Xp  Z  _ 

■Az 

0 

The  characteristic  equations  results  from  setting  the  determinant  of  the 
.coefficient  matrix  to  zero 


(k-  -  I  CO^)  (k,  -  I  0)2)  -  O)^  =  0 

6  y  ^  yz 

(I  I  -  l2  Jo/  -  (k.l  +  k,  I  )o;2  +  k  k,  =  0 
y  z  yZ  6  z  ij/y  6  ^ 

let  l2  =  I  I  - 

yzyz 

then,  the  natural  frequencies  of  the  system  are 


(A-7) 


(A- 8) 


0) 


xly,6 


k  I  +  k  I 
6  z  u/  y 
2I2 


+ 


k  I  +  k  I  I 

9  ^  ^  v\ 

2i2  ! 


k-k, 

e  \j/ 

l2 


(A-9) 


The  amplitude  ratio  for  the  frequency  O) 


0)^  into  Eq.  (A-7) 


is  obtained  by  substituting 


-  I  <^2 

9  y  yp 

I  0)2  1 

yz  Xp 

j 

_  1 

0 

I  0)? 

L  yz  i// 

k  -  I  0)2 
xy  z  i//J 

1 

0 

(k.  -  I  0)f  )Ai  ,  +  I  0)2  A,.  =  0 

0  y  »//  V  yz  xf^ 


and, 


yz 


xj/  V 


- 


I  w2)A2, 

Z  x}/ 


then. 


Ai  ,  I  0)^. 

^  _  yz  i// 

^2^  %'\^X() 


(A-10) 


The  second  subscript  on  A  is  used  to  indicate  that  the  amplitudes  are 


associated  with  the  frequency  0)^.  Thuj},|One  possible  solution  for  the 
equations  of  motion  is  given  by  the  trial  solution,  Eq.  (A-6),  in  which  the 
frequency  is  and  the  amplitudes  are  related  by  Eq.  (A-10).  Another 

solution  is  possible  if  sin  cot  is  replaced  by  cos  cot,  i.  e.  , 


cos  cot 


A  procedure  similar  to  that  for  Eq.  (A-^  gives  frequencies  identical  to  that 


9 

1  1 

Bi 

r  i 

Bz 

(A-11) 


of  Eq.  (A- 9).  and  an  amplitude  ratio  for  co^  as 


Bj  ,  A 1  ,  I  CO^. 

- 

Bt.,  At.,  k.-I  CO; 

^  ^  d  y  x}y 


(A-12) 


The  complete  solution  for  a  free  vibration  at  the  first  natural  frequency, 


CO^  ,  is  given  by 


6  ~  Ai  ,  sinco,  t  +  Bi  ,  cos  CO,  t 

P  ^  \ly 

\1/  =  At,  sinco,  t  +  Bt,  cos  CO,  t 
^  \jy  \}/ 


(A-13) 

(A-14) 


A  similar  procedure  for  the  second  natural  frequency,  CO,,  ,  yields 

V 


6  =  Ai  sinco  t  +  Bi  cos  CO  t 
U .  U  U  V 

^  =  Az^  sinCO^t  +  B2^  cos  CO^t 


(A-15) 

(A-16) 


in  which 


Zle. 


B 


Be. 


I  col  k,-  I  CO^ 

yz  6  _  xy  z  6 

k  -  I  CO^  ~  I  CO^ 

9  y  9  yz  9 


(A-17) 
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The  most  general  solution  for  the  free  vibrations  of  the  system  is  given 


by  the  superposition  of  the  motions  in  the  two  principal  modes  of  vibration, 
given  by  Eqs.  (A-13)  to  (-A-I6). 


6  =  Ai  ,  sinCtJ.t  +  Bi  ,  cosCU.t  +  Ai  ^sincU-t  +  Bi-cosco.t 

ij/  xl/  l//  xjy  6  6  9  6 

^  +-A2gsinw^t  +  Bz^cosco^t 


(A-18) 

(A-19) 


'making  use  of  Eqs.  (A-12)  arid  (A-17) 


1  r 

_ hk 


\  A —  +  B2aCOsW^t^(A-20) 


^9 - 9  "9 - 6 


k,,  sinCU,  t  +  Bz  ,  cosCO,  t  >  +  ■(Az.^sinai^t  +  B2„cosa),,t 

V  ^  V  \k  \  )  0  0  0  & 


(A-21) 


Eqs.  (A-20)  and  (A-21)  can  be  written  as 


0  =  .011  qi  +  0i2qE 


(A-22) 


\k  =  921  qi  +  ^2^2 


(A-23) 


where, 


I  w; 

vz  1/ 


I  017  -  k, 
Z  XU  XL 


11  I  CO,  -k^  I  co; 
y  xj/  6  yz  0/ 


(A -24) 


I  co;  I  co^^'-k, 

A  yz  9  _  z  9  0 

*^12  “  I  co-^-k  “  I  co^ 
y  6  6  yz  9 


(A-25) 


02J  =  922  =  1 


(A-26) 


Qi  =  Az,  sinco.t  +  Bz,  cosCO  t 
0/0/  0/0/ 


{A-27) 


q2=  Az^  sinco^t  +  Bz^  cosco^t 


(A-28) 


The  constants  Az,  »  Bz,  ,  Az  ,  and  Bz^  can  be  determined  from  the 

V  0/  0  0 


initial  conditions 


(A-29) 


0(0)  =  e 


0(0)  =  0 


ll^iO)  =  »//(0)=  1//^ 


(A-30) 


The  quantities  qj  and  qz  represent  the  amplitudes  of  the  motion  in 
the  two  normal  modes  and  are  called  the  normal  coordinates.  A  work-and 
-energy  approach  to  the  equations  of  motion  shows  that  the  multiplication  of 
the  equations  of  motion  in  the  normal  coordinate  system  by  the  constant 
will  uncouple  the  equations  of  motion. 

Now  consider  the  complete  equations  of  motion,  Eqs.  (A-1)  and  (A-2), 
which  can  be  written  as  . 


l9l2  ij 


■  I  -  I  “ 

0 

0  1 

0 

(K  o' 

0 

0 

y  yz 

•  • 

+ 

0 

+ 

0 

~ 

-I  I 

_  yz  z  _ 

o 

_ 1 

0^ 

0  k, 

L  i^J 

1 

0 

(A-31) 


From  Eqs.  (A-22)  and  (A-23) 


0 

r0 

^11 

0 

^12 

qi 

-1 

1  - 

q2 

Substituting  Eq.  (A- 32)  into  Eq.  (A-31)  and  multiplying  by 

uncoupled  equations  in  the  normal  coordinate  system 

“I 


1' 

[912  1 


(7)  1 

ai 

(1)  •  1 
az- 


y  yz 

I  I 
L  yz  z  , 


0  0 

11  12 

1  1 


qi 

qz 


+ 


0  1 

^11 


0  - 


^0  “ 


0  c 


will  give 


011  012 
1  1 


0, 


11 


0, 


12 


^0  ^ 


0  k 


x[y. 


0  0 

11  12 

1  1 


(A-32) 


matrix  multiplication  gives 


^7 


le  +k, 

,  ^11  e  xjy 

'  77^ 

CO, 


(I)  ^  k.  +k , 
^12  0^  xp 

0)2 

e 


qi 


qa) 


S  +  S/. 


f  C„+C, 

12  9  ^ 


qi 


\  qz 


qi 

1 

= 

0 

qz 

0 

''■..“e'V  ° 

“  >.2‘'s+V 

where  it  was  assumed  tha"  ((h  0,^C^  +  C,  )  is  much  smaller  than  either 

'^11  ^iz  Q  y/  ' 

^^12^9  S//^’  ^9  ^xjj  givea  by  Eq.  (A-9) 

Equation  (A-33)is  uncoupled  and  can  be  rewritten  as 


qi  + 


qz  + 


02  k  +k, 
qz  0  xL- 


^q^2  +  i^^qz  =  0 


Equations  (A-34)  and  (A-35)  can  be  solved  in  the  usual  manner  to  give 

-  61//CO.  1 1 


qi=  e  W 
-a«  co^t 


qz  =  e 


’Q 


A,  cosCO,  t  +  B,  sinco,  t 

^  ^  i' 


A  coscO  t  +  sincO-t 
9  9  9  9 


where 


and, 


V 


^  c„  +  C, 
q  1  9  xp 

f  k.+  k, 
qi  9  xjj 


CO, 

J' 


9  ^n+  C, 

qz  0  xjf 


0)^  k^  +  k, 

qa  0  l// 


CO 


0 


The  results,  Eqs.  (A-36)  and  (A-37)  can  now  be  transformed  into  the  0  and  l// 
coordinates  using  Equation  (A-32). 


(A-33) 


(A-34) 


(A-35) 


(A-36) 

(A-37) 

(A-38) 


(A-39) 


28 


0(t)  =  <!>,,  ^ 

\{/(t)  =  |A^cosa)^t+  B^sinCO^tje 
+  lAgCosco^t  +  Bg  sincu^t  le  ^6^ 


X.i.t 


(A-40) 


(A-41) 


where  X,=  6,0),  and  X^  =  6^0). 

;//;//  i//  9  9  6 

The  constants  in  Eqs.  (A-40)  and(A-41)  can  be  determined  from  the  initial 
conditions,  Eqs.  (A-29)  and  (A-30). 


(A-42) 


A,  = 


6  -  (i> 

O  >2  o 


6  -  <f>  4/ 

B,  =  -V^-—  +  5,  A, 

4/  ^O)  4y  4/ 


(b  4^-9 

A  -  u  ° _ ° 


■^9  = 


(ji  4^  -  9 

-°-  +  6-A 

(j)(jO^  9  9 


where  6=6  -6 

1  ^11  h2 


(A -4  3) 


The  solution  for  the  damped  motion,  coupled  by  is 


6(t)  =  ^  e 

12 


+  e 

11 


■  V  jft.  %-  “c 


>  .  (  /  0  -(/>  4> 
3  (  o  /u  o 


cosa).t  + 

C7 


+ 


^  /{fi  if/  -o  \ 


(^4 


'0  ^9 


!> 


COSO),  t  + 

4^ 


r9  -6  'P  X,  ,6  -  6  i/^\' 
I  O  ^2^0  ^  (  O  ,z^o\ 


!>  ) 


sino^tj 


s  ino) ,  t 
4  I 


(A-44) 


.\/9  -6  4^ 

*  \[  ^ 


_L2_ 


V-0) 


+ 


0)., 


_^2_£ 

sino),  t 

\  /  _ 

+  e 


0 


cosw  t  + 
u 


4,  > 

•  « 

'6  4^  -9  X  /6  4>  -9 
1 1  o  o  ,  _0.  f  11  O  o 


sinOlgt 


(A-45) 


where  0,0  and  are  given  by  Eqs.  (A-24),  (A-25)  and  (A-43)  respectively. 
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